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This paper gives a dichotomy theorem for the complexity of computing the partition 
function of an instance of a weighted Boolean constraint satisfaction problem. The prob- 
lem is paramctcrised by a finite set T of non-negative functions that may be used to 
assign weights to the configurations (feasible solutions) of a problem instance. Classical 
^sO ■ constraint satisfaction problems correspond to the special case of 0,1-valued functions. 

We show that computing the partition function, i.e. the sum of the weights of all config- 
■^J- ■ urations, is FP* P -complete unless either (1) every function in T is of "product type", or 

(2) every function in T is "pure affine" . In the remaining cases, computing the partition 
function is in P. 

O 

> : 

■ 1 Introduction 

This paper gives a dichotomy theorem for the complexity of the partition function of weighted 
Boolean constraint satisfaction problems. Such problems are parameterised by a set J- of 
non-negative functions that may be used to assign weights to configurations (solutions) of 
the instance. These functions take the place of the allowed constraint relations in classical 
constraint satisfaction problems (CSPs). Indeed, the classical setting may be recovered by 
restricting T to functions with range {0, 1}. The key problem associated with an instance 
of a weighted CSP is to compute its partition function, i.e., the sum of weights of all its 
configurations. Computing the partition function of a weighted CSP may be viewed a gen- 
eralisation of counting the number of satisfying solutions of a classical CSP. Many partition 
functions from statistical physics may be expressed as weighted CSPs. For example, the 
Potts model [23] is naturally expressible as a weighted CSP, whereas in the classical frame- 
work only the "hard core" versions may be directly expressed. (The hard-core version of the 
antiferromagnetic Potts model corresponds to graph colouring and the hard-core version of 
the ferromagnetic Potts model is trivial — acceptable configurations colour the entire graph 
with a single colour.) A corresponding weighted version of the decision CSP was investigated 
by Cohen, Cooper, Jeavons and Krokhin |3' . This results in optimisation problems. 
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We use #CSP(JF) to denote the problem of computing the partition function of weighted 
CSP instances that can be expressed using only functions from T . We show in Theorem [J] 
below that if every function / € T is "of product type" then computing the partition function 
Z(I) of an instance I can be done in polynomial time. Formal definitions are given later, 
but the condition of being "of product type" is easily checked — it essentially means that 
the partition function factors. We show further in Theorem 0] that if every function / € T is 
"pure affine" then the partition function of Z(I) can be computed in polynomial time. Once 
again, there is an algorithm to check whether T is pure affine. For each other set T , we show 
in Theorem H] that computing the partition function of a #CSP(JF) instance is complete for 
the class FP* P . The existence of algorithms for testing the properties of being purely affine 
or of product type means that the dichotomy is effectively decidable. 

1.1 Constraint satisfaction 

Constraint Satisfaction, which originated in Artificial Intelligence, provides a general frame- 
work for modelling decision problems, and has many practical applications. (See, for exam- 
ple [18J.) Decisions are modelled by variables, which are subject to constraints, modelling 
logical and resource restrictions. The paradigm is sufficiently broad that many interesting 
problems can be modelled, from satisfiability problems to scheduling problems and graph- 
theory problems. Understanding the complexity of constraint satisfaction problems has be- 
come a major and active area within computational complexity [71 114j. 

A Constraint Satisfaction Problem (CSP) typically has a finite domain, which we will 
denote by [q] = {0, 1 . . . , q — 1} for a positive integer q^ A constraint language T with 
domain [q] is a set of relations on [q]. For example, take q = 2. The relation R = {(0,0, 1), 
(0, 1,0), (1, 0, 0), (1, 1, 1)} is a 3-ary relation on the domain {0, 1}, with four tuples. 

Once we have fixed a constraint language T, an instance of the CSP is a set of variables 
V = {vi, . . . ,v n } and a set of constraints. Each constraint has a scope, which is a tuple of 
variables (for example, (v^,v^,vi)) and a relation from T of the same arity, which constrains 
the variables in the scope. A configuration a is a function from V to [q]. The configuration a 
is satisfying if the scope of every constraint is mapped to a tuple that is in the corresponding 
relation. In our example above, a configuration a satisfies the constraint with scope (t>4, v§, v\) 
and relation R if and only if it maps an odd number of the variables in {fi,f4,fs} to the 
value 1. Given an instance of a CSP with constraint language V, the decision problem CSP(r) 
asks us to determine whether any configuration is satisfying. The counting problem #CSP(r) 
asks us to determine the number of (distinct) satisfying configurations. 

Varying the constraint language V defines the classes CSP and #CSP of decision and 
counting problems. These contain problems of different computational complexities. For 
example, if T = {Ri, R2, R3} where R±, R2 and R3 are the three binary relations defined 
by R x = {(0,1), (1,0), (1,1)}, R 2 = {(0,0), (0,1), (1,1)} and R 3 = {(0, 0), (0, 1), (1, 0)}, then 
CSP(r) is the classical 2-Satisfiability problem, which is in P. On the other hand, there 
is a similar constraint language V with four relations of arity 3 such that 3-Satisfiability 
(which is NP-complete) can be represented in CSP(r'). It may happen that the counting 
problem is harder than the decision problem. If T is the constraint language of 2-Satisfiability 
above, then #CSP(r) contains the problem of counting independent sets in graph, and is 
#P-complete [22], even if restricted to 3-regular graphs [12]. 

1 Usually [q] is defined to be {1, 2, ... , q}, but it is more convenient here to start the enumeration of domain 
elements at rather than 1. 
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Any decision problem CSP(r) is in NP, but not every problem in NP can be represented 
as a CSP. For example, the question "Is G Hamiltonian?" cannot naturally be expressed 
as a CSP, because the property of being Hamiltonian cannot be captured by relations of 
bounded size. This limitation of the class CSP has an important advantage. If P 7^ NP, 
then there are problems which are neither in P nor NP-complete [TB]. But, for well-behaved 
smaller classes of decision problems, the situation can be simpler. We may have a dichotomy 
theorem, partitioning all problems in the class into those which are in P and those which 
are NP-complete. There are no "leftover" problems of intermediate complexity. It has been 
conjectured that there is a dichotomy theorem for CSP. The conjecture is that CSP(r) is 
in P for some constraint languages T, and CSP(T) is NP-complete for all other constraint 
languages T. This conjecture appeared in a seminal paper of Feder and Vardi [10], but has 
not yet been proved. 

A similar dichotomy, between FP and #P-complete, is conjectured for f^CSP [2]. The 
complexity classes FP and ffP are the analogues of P and NP for counting problems. FP 
is simply the class of functions computable in deterministic polynomial time. #P is the 
class of integer functions that can be expressed as the number of accepting computations of 
a polynomial-time non-deterministic Turing machine. Completeness in #P is defined with 
respect to polynomial-time Turing reducibility [T71 Chap. 18]. Bulatov and Dalmau [2] have 
shown in one direction that, if f^CSP(r) is solvable in polynomial time, then the constraints 
in r must have certain algebraic properties (assuming P 7^ #P). In particular, they must 
have a so-called Mal'tsev polymorphism. The converse is known to be false, though it remains 
possible that the dichotomy (if it exists) does have an algebraic characterisation. 

The conjectured dichotomies for CSP and #CSP are major open problems for computa- 
tional complexity theory. There have been many important results for subclasses of CSP and 
#CSP. We mention the most relevant to our paper here. The first decision dichotomy was 
that of Schaefer [19] . for the Boolean domain {0, 1}. Schaefer's result is as follows. 

Theorem 1 (Schaefer pj|]). Let T be a constraint language with domain {0, 1}. The problem 
CSP(r) is in P ifT satisfies one of the conditions below. Otherwise, CSP(r) is NP-complete. 

(1) r is 0-valid or 1-valid. 

(2) r is weakly positive or weakly negative. 

(3) r is affine. 

(4) r is bijunctive. 

We will not give detailed definitions of the conditions in Theorem [TJ but the interested 
reader is referred to the paper [19] or to Theorem 6.2 of the textbook [7] . An interesting feature 
is that the conditions in [7J Theorem 6.2] are all checkable. That is, there is an algorithm 
to determine whether CSP(r) is in P or NP-complete, given a constraint language T with 
domain {0, 1}. Creignou and Hermann [6] adapted Schaefer's decision dichotomy to obtain a 
counting dichotomy for the Boolean domain. Their result is as follows. 

Theorem 2 (Creignou and Hermann [6]). Let T be a constraint language with domain {0, 1}. 
The problem #CSP(r) is in FP if T is affine. Otherwise, #CSP(r) is ffV-complete. 

A constraint language T with domain {0, 1} is affine if every relation R E V is affine. A 
relation R is affine if the set of tuples x E R is the set of solutions to a system of linear 
equations over GF(2). These equations are of the form v\ © • • • © v n = and V\ © • • • © v n = 1 
where © is the exclusive or operator. It is well known (see, for example, Lemma 4.10 of [7J) 
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that a relation R is affine iff a, b, c £ R implies d = a©6©c£i?. (We will use this 
characterisation below.) There is an algorithm for determining whether a Boolean constraint 
language T is affine, so there is an algorithm for determining whether #CSP(r) is in FP or 
#P-complete. 

1.2 Weighted #CSP 

The weighted framework of [1] extends naturally to Constraint Satisfaction Problems. Fix 
the domain [q]. Instead of constraining a length-A; scope with an arity-fe relation on [q], 
we give a weight to the configuration on this scope by applying a function / from [q] k to 
the non-negative rationals. Let T q = {/ : [q] k — ► Q + | G N} be the set of all such 
functions (of all arities)H Given a function / G J- q of arity /c, the underlying relation of / 
is given by ify = {x G [g] fc | f(x) / 0}. It is often helpful to think of R f as a table, with 
k columns corresponding to the positions of a £;-tuple. Each row corresponds to a tuple 
x = (xi, . . . , Xk) G Rf. The entry in row x and column j is Xj, which is a value in [q]. 

A weighted #CSP problem is parameterised by a finite subset J 7 of T q , and will be 
denoted by #CSP(J r ). An instance / of #CSP(J r ) consists of a set V of variables and a 
set C of constraints. Each constraint C G C consists of a function /c G .F (say of arity A:c) 
and a scope, which is a sequence sc = {vc,i, ■ ■ ■ ,vc,k c ) °f variables from V. The variables 
vc,i, ■ ■ ■ , vc,k c n eed not be distinct. As in the unweighted case, a configuration a for the 
instance I is a function from V to [q]. The weight of the configuration a is given by 

w{cr) = Yl fc(o-(v c ,i), (T{v c ,k c ))- 
Cec 

Finally, the partition function Z(I) is given, for instance /, by 

z(i)= y, a) 

<rV-*[q] 

In the computational problem #CSP(J r ), the goal is to compute Z(I), given an instance /. 

Note that an (unweighted) CSP counting problem #CSP(r) can be represented naturally 
as a weighted CSP counting problem. For each relation R £ T, let f R be the indicator 
function for membership in R. That is, if x G R we set f R {x) = 1. Otherwise we set 
f R {x) = 0. Let J- = {f R | R G r}. Then for any instance / of #CSP(r), the number of 
satisfying configurations for / is given by the (weighted) partition function Z{I) from (JTJ) . 

This framework has been employed previously in connection with graph homomorphisms pQ. 
Suppose H = (Hij) is any symmetric square matrix H of rational numbers. We view H as 
being an edge-weighting of an undirected graph 7i, where a zero weight in H means that 
the corresponding edge is absent from 7i. Given a (simple) graph G = (V, E) we consider 
computing the partition function 

Zh{G)= Y u '( (T )' wn ere w(a) = \\ H <r(u)a(v)- 

a:V^[q] {u,v}eE 

Within our framework above, we view H as the binary function h : [q] 2 — ► R, and the problem 
is then computing the partition function of #CSP({h}). 

2 We assume £ N, so we allow non-negative constants. 



4 



Bulatov and Grohe [I] call H connected if Ti is connected and bipartite if TL is bipartite. 
They give the following dichotomy theorem for non-neg ative HE 

Theorem 3 (Bulatov and Grohe [4]). Let H be a symmetric matrix with non-negative rational 
entries. 

(1) If H is connected and not bipartite, then computing Zh is in FP if the rank of H is at 
most 1; otherwise computing Zh is ffP-hard. 

(2) If H is connected and bipartite, then computing Zu is in FP if the rank of H is at most 2; 
otherwise computing Zh is #P-hard. 

(3) If H is not connected, then computing Zh is in FP if each of its connected components 
satisfies the corresponding conditions stated in (1) or (2); otherwise computing Zh is 
#P-hard. □ 

Many partition functions arising in statistical physics may be viewed as weighted #CSP 
problems. An example is the q-state Potts model (which is, in fact, a weighted graph ho- 
momorphism problem). In general, weighted #CSP is very closely related to the problem of 
computing the partition function of a Gibbs measure in the framework of Dobrushin, Lanford 
and Ruelle (see [T] ) . See also the framework of Scott and Sorkin [20J . 

1.3 Some Notation 

We will call the class of (rational) weighted #CSP problems weighted #CSP. The sub-class 
having domain size q = 2 will be called weighted Boolean #CSP, and will be the main focus 
of this paper. We will give a dichotomy theorem for weighted Boolean #CSP. 

Since weights can be arbitrary non-negative rational numbers, the solution to these prob- 
lems is not an integer in general. Therefore #CSP(JF) is not necessarily in the class #P. 
However, Goldberg and Jerrum [11] have observed that Z(I) = Z(I)/K(I), where Z is a 
function in #P and K(I) is a positive integer computable in FP. This follows because, for all 
/ £ T , we can ensure that /(•) = f(-)/K(I), where /(•) 6 N, by "clearing denominators". The 
denominator K(I) can obviously be computed in polynomial time, and it is straightforward 
to show that computing Z{I) is in #P, so the characterisation of [TT] follows. The resulting 
complexity class, comprising functions which are a function in #P divided by a function in 
FP, is named #Pq in [TT], where it is used in the context of approximate counting. Clearly 
we have 

weighted #CSP C #P Q C FP #P . 

On the other hand, if Z(I) £ weighted #CSP is ^P-hard, then, using an oracle for computing 
Z(I), we can construct a #P oracle Z(I) as outlined above. (Note that Z{I) ^ #P in 
general.) Using this, we can compute any function in FP* P with a polynomial time-bounded 
oracle Turing machine. Thus any #P-hard function in weighted #CSP is complete for FP* P . 
We will use this observation to state our main result in terms of completeness for the class 
FP# P . 

We make the following definition, which relates to the discussion above. We will say that 
F Q Fq simulates f € T q if, for each instance I of #CSP(J 7 U {/}), there is a polynomial time 
computable instance V of jfdSP[T\ such that Z(I) = tp(I)Z(I') for some cp(I) € Q which is 
FP-computable. This generalises the notion of parsimonious reduction [J7| among problems in 

3 This is not quite the original statement of the theorem. We have chosen here to restrict all inputs to be 
rational, in order to avoid issues of how to represent, and compute with, arbitrary real numbers. 
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ffP. We will use <t to denote the relation "is polynomial-time Turing-reducible to" between 
computational problems. Clearly, if T simulates /, we have #CSP(J r U {/}) <t #CSP(JT). 
Note also that, if / = Kf, for some constant K > 0, then {/} simulates /. Thus there is no 
need to distinguish between "proportional" functions. 

We use the following terminology for certain functions. Let x= be the binary equality 
function defined on [q] as follows. For any element c G [q], x=( c i c) = 1 and for any pair (c, d) 
of distinct elements of [q], x=( c id) = 0. Let \=£ be the binary disequality function given by 
X^(c, <£) = ! — x=( c > d) for all c, d G [g]o We say that a function / is of product type if / can 
be expressed as a product of unary functions and binary functions of the form x= an d x^- 

We focus attention in this paper on the Boolean case, q = 2. In this case, we say that a 
function / G Ti has affine support if its underlying relation Rf, defined earlier, is affine. We 
say that / is pure affine if it has affine support and range {0, w} for some w > 0. Thus a 
function is pure affine if and only if it is a positive real multiple of some (0,1-valued) function 
which is affine over GF(2). 

1.4 Our Result 

Our main result is the following. 

Theorem 4. Suppose T " C T 2 = {/ : {0, l} fe -> Q+ | k G N}. If every function in T is of 
product type then f£CSP(J-) is in FP. If every function in J- is pure affine then #CSP(J-) is 
in FP. Otherwise, #CSP(.F) is FP* P -complete. 

Proof. Suppose first that J- is of product type. In this case the partition function Z{I) of 
an instance I with variable set V is easy to evaluate because it can be factored into easy- 
to-evaluate pieces: Partition the variables in V into equivalence classes according to whether 
or not they are related by an equality or disequality function. (The equivalence relation on 
variables here is "depends linearly on".) An equivalence class consists of two (possibly empty) 
sets of variables U\ and C/ 2 . All of the variables in U\ must be assigned the same value by a 
configuration a of nonzero weight, and all variables in U2 must be assigned the other value. 
Variables in U1UU2 are not related by equality or disequality to variables inV\(£/iU£/2)- The 
equivalence class contributes one weight, say a, to the partition function if variables in U\ 
are given value "0" by a and it contributes another weight, say (3, to the partition function 
if variables in U\ are given value "1" by a. Thus, Z(T) = {a + f3)Z(I'), where I' is the 
instance formed from / by removing this equivalence class. Therefore, suppose we choose any 
equivalence class and remove its variables. Since J- contains only unary, equality or binary 
disequality constraints, we can also remove all functions involving variables in U\ U XJ% to give 
J-'. Then /' is of product type with fewer variables, so we may compute Z{I') recursively. 

Suppose second that T if pure affine. Then Z(T) = Ylf^w^ Z(I'), where {0,wj} is the 
range of /, kf is the number of constraints involving / in I, and I' is the instance obtained 
from / by replacing every function / by its underlying relation Rf (viewed as a function with 
range {0,1}). Z(I') is easy to evaluate, because this is just counting solutions to a linear 
system over GF(2), as Creignou and Hermann have observed [6]. 

Finally, the #P-hardness in Theorem [4] follows from Lemma [5] below. □ 

Lemma 5. If f G Ti is not of product type and g(zJ~2 is not pure affine then #CSP({/, g}) 
is ffP-hard. 

4 A more general disequality function is defined in the Appendix. 
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Note that the functions / and g in Lemma [5] may be one and the same function. So 
#CSP({/}) is #P-hard when / is not of product type nor pure affine. The rest of this article 
gives the proof of Lemma [SJ 

2 Useful tools for proving hardness of #CSP 

2.1 Notation 

For any sequence ui, . . . ,uj~ of variables of I and any sequence ci, . . . ,c/fe of elements of the 
domain [q], we will let Z(I \ o~(u\) = ci, . . . , o~(uk) = Ck) denote the contribution to Z(I) from 
assignments a with a(u\) = c±, ■ ■ ■ , a(uk) = c&. 

2.2 Projection 

The first tool that we study is projection, which is referred to as "integrating out" in the 
statistical physics literature. 

Let / be a function of arity k, and let J = . . . ,j r } be a size-r subset of {1, . . . , k}, 
where ji < • • • < j r u We say that a /c-tuple x' € [q] k extends an r-tuple x € [q] r on J (written 
x' 3 J x) if x' agrees with x on indices in J; that is to say, x'-. = X{ for all 1 < i < r. The 
projection g of / onto J is defined as follows. For every x G [q] r , g(x) = Ylx'z\jx f( x ')- 

The following lemma may be viewed as a weighted version of Proposition 2 of [2], where 
it is proved for the unweighted case. It is expressed somewhat differently in [2], in terms of 
counting the number of solutions to an existential formula. 

Lemma 6. Suppose J- C T q . Let g be a projection of a function f £ J- onto a subset of its 
indices. Then #CSP(J r U {g}) < T #CSP(J r ). 

Proof. Let k be the arity of / and let g be the projection of / onto the subset J of its indices. 
Let I be an instance of #CSP(JFU {g}). We will construct an instance I' of #CSP(^-") such 
that Z(I) = Z(I'). The instance I' is identical to I except that every constraint C of I 
involving g is replaced with a new constraint C of I 1 involving /. The corresponding scope 
(yc',1, • • ■ , vc,k) is constructed as follows. If ji is the th element of J, then v' c , ^ = vc : £- The 
other variables, vcj (j ^ J), are distinct new variables. We have shown that T simulates g 
with (f)(1) = 1. □ 

2.3 Pinning 

For c G [q], 5 C denotes the unary function with <5 c (c) = 1 and S c (d) = for d ^ c. The 
following lemma, which allows "pinning" CSP variables to specific values in hardness proofs, 
generalises Theorem 8 of [2], which does the unweighted case. Again [2] employs different 
terminology, and its theorem is a statement about the full idempotent reduct of a finite 
algebra. The idea of pinning was used previously by Bulatov and Grohe of [1] in the context 
of counting weighted graph homomorphisms (see Lemma 32 of [4]). A similar idea was used 
by Dyer and Greenhill in the context of counting unweighted graph homomorphisms — in 
that context, Theorem 4.1 of |8j allows pinning all variables to a particular component of the 
target graph H. 

5 It is not necessary to choose this particular ordering for J, but it is convenient to do so. 
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Lemma 7. For every T C T q , #CSP(F U \J ce[q] S c ) < T #CSP(.F). 

The proof of Lemma [7] is deferred to the appendix. Since we only use the case q = 2 in 
this paper, we provide the (simpler) proof for the Boolean case here. 

Lemma 8. For every T C T 2 , #CSP(T U {S , <5i}) < T #CSP(.T). 

Proof. For a; G [2] fc , let x be the fc-tuple whose i'th component, Xj, is Xj © 1, for all i. Say 
that T is symmetric if it is the case that for every arity-A; function / G T and every x £ [2] fc , 

m = f(x). 

Given an instance / of #CSP(J r U {<5o, #i}) with variable set V we consider two instances 
/' and I" of ^CSP(JT). Let Vo be the set of variables v of I to which the constraint 5o(v) is 
applied. Let V\ be the set of variables v of I to which the constraint 5i(v) is applied. We 
can assume without loss of generality that Vo and V\ do not intersect. (Otherwise, Z(I) = 
and we can determine this without using an oracle for #CSP(JT).) Let Vi = V \ (Vo U VI). 
The instance /' has variables V2 U {to,ti} where to an d t\ are distinct new variables that are 
not in V. Every constraint C of / involving a function / 6 T corresponds to a constraint C 
of C is the same as C except that variables in Vo are replaced with to and variables 
in Vi are replaced with t\. Similarly, the instance I" has variables V2 U {t} where t is a new 
variable that is not in V. Every constraint C of / involving a function / E T corresponds to 
a constraint C" of I". The constraint C" is the same as C except that variables in Vo U Vi 
are replaced with t. 

Case 1. J- is symmetric: By construction, 

Z(I') - Z(I") = Z(I' I <r(io) = 0,<r(ti) = 1) + Z(I' I ff (t ) = l )( r(ti) = 0). 
By symmetry, the summands are the same, so 

Z(I') - Z(I") = 2Z(I' I a(t ) = 0,a(h) = 1) = 2Z{I). 

Case 2. T is not symmetric: Let / be an arity-/c function in T and let x G [2] k so that 
/(x) > /(x) > 0. Let s = (t xi , . . . ,t Xk ) and let I' x be the instance derived from I' by adding a 
new constraint with function / and scope s. Similarly, let I" be the instance derived from /" 
by adding a new constraint with function / and scope (t, . . . ,t). Now 

Z(I' X ) = Z(I' I a(t ) = 0,a(h) = \)f{x) + Z{l' \ a(t ) = l,a(h) = 0)/(x) 

+ Z(l' I a(t ) = 0, a(h) = 0)/(0, . . . , 0) + Z(I' \ a(t ) = 1, a(h) = 1)/(1, . . . , 1) 
= Z(I' I a(t ) = 0, a(h) = l)/(x) + Z(I' \ a(t ) = 1, a(h) = 0)/(x) + Z(/*). 

Thus we have two independent equations, 

Z(4) - Z(/*) = Z(/' I a(t ) = 0, aCtj) = l)/(x) + Z(I' I a(t ) = 1, a^) = 0)/(x), 
Z(I')-Z(I") = Z(I' \a{to) = 0,a(h) = l) + Z(I' \ a(t ) = l,a(h) = 0) 

in the unknowns Z{1' \ cr(to) = 0,cr(ti) = 1) and Z{I' \ a(to) = IjC^l) = 0). Solving these, 
we obtain the value of Z(I' \ a(to) = 0, cr(ti) = 1) = Z(I). □ 
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2.4 #P-hard problems 



To prove Lemma [SJ we will give reductions from some known #P-hard problems. The first of 
these is the problem of counting homomorphisms from simple graphs to 2-vertex multigraphs. 
We use the following special case of Bulatov and Grohe's Theorem [3l 

Corollary 9 (Bulatov and Grohe [4]). Let H be a symmetric 2x2 matrix with non-negative 
real entries. If H has rank 2 and at most one entry of H is then Eval(H) is #P-hard. 

We will also use the problem of computing the weight enumerator of a linear code. Given a 
generating matrix A G {0, l} rxC ' of rank r, a code word c is any vector in the linear subspace T 
generated by the rows of A over GF(2). For any real number A, the weight enumerator of the 
code is given by Wa(A) = 2~^ ce x^" C "> wnere ll c ll is the number of l's in c. The problem of 
computing the weight enumerator of a linear code is in FP for A G {— 1,0,1}, and is known to 
be #P-hard for every other fixed A G Q (see [23]). We could not find a proof, so we provide 
one here. We restrict attention to positive A, since that is adequate for our purposes. 

Lemma 10. Computing the Weight Enumerator of a Linear Code is fj^P-hard for any fixed 
positive rational number A ^ 1. 

Proof. We will prove hardness by reduction from a problem EVAL(-ff), for some appropriate H, 
using Corollary [9j Let the input to EvAL^iJ) be a connected graph G = (V,E) with V = 
{vi, . . . , v n } and E = {e\, . . . , e m }. Let B be the n x m incidence matrix of G, with bij = 1 
if Vi G €j and bij = otherwise. Let A be the (n — 1) x m matrix which is B with the row 
for v n deleted. A will be the generating matrix of the Weight Enumerator instance, with 
r = n — 1 and C = m. It has rank (n — 1) since G contains a spanning tree. A code word 
c has Cj = @i e jjbij, where U C V \ {v n }. Thus Cj = 1 if and only if ej has exactly one 
endpoint in U, and the weight of c is A fc , where k is the number of edges in the cut U,V\U. 
Thus W^(A) = ^Zjj(G), where H is the symmetric weight matrix with Hn = H22 = 1 and 
H\2 = H21 = A. The 5 arises because we fixed which side of the cut contains v n . Now H 
has rank 2 unless A = 1, so this problem is #P-hard by Corollary [9j Note, by the way, that 
Zh{G) is the partition function of the Ising model in statistical physics [5]. □ 

3 The Proof of Lemma [5] 

Throughout this section, we assume q = 2. The following Lemma is a generalisation of a result 
of Creignou and Hermann [6], which deals with the case in which / is a relation (or, in our 
setting, a function with range {0, 1}). The inductive technique used in the proof of Lemma [TT1 
(combined with the follow-up in Lemma [T2l) is good for showing that jfCSP(J-) is #P-hard 
when T contains a single function. A very different situation arises when #CSP({/}) and 
#CSP({<?}) are in FP but #CSP({/, g}) is #P-hard due to interactions between / and g - 
we deal with that problem later. 

Lemma 11. Suppose that f G J-% does not have affine support. Then #CSP({/}) is jf-P-hard. 

Proof. Let k be the arity of /, and let us denote the i th component of fc-tuple a G Rf by a«. 
The proof is by induction on k. The lemma is trivially true for k = 1, since all functions of 
arity 1 have affine support. 
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For k = 2, we note that since Rf is not affine, it is of the form Rf = {(a, (3), (at, f3), (at, (3)} 
for some a € {0, 1} and (3 £ {0, 1}. We can show that #CSP({/}) is #P-hard by reduction 
from Eval(H) using 

/ /(0,0) /(0,1) \ 

V /(i.o) f(i,i) )' 

which has rank 2 and exactly one entry that is 0. Given an instance G = (V,E) of Eval(H) 
we construct an instance / of #CSP({/}) as follows. The variables of / are the vertices of G. 
For each edge e = (u, v) of G, add a constraint with function / and variable sequence u, v. 
Corollary H now tells us that EvAh(H) is #P-hard, so #CSP({/}) is #P-hard. 

Suppose k > 2. We start with some general arguments and notation. For any i £ 
{1, . . . ,k} and any a £ {0, 1} let f t=a be the function of arity k — 1 derived from / by pinning 
the i'th position to a. That is, f l=a (x\, . . . , x^-i) = f(x\, . . . , Xi-%, a, Xi+\, . . . , x^). Also, 
let f t= * be the projection of / onto all positions apart from position i (see Section 12 . 2[) . 
Note that #CSP({/ i=a }) < T #CSP({/, <5 , 5{\), since f i=a can obviously be simulated by 
{/,<Jo,<5i}. Furthermore, by Lemma #CSP({/, 6 , 5i}) < T #CSP({/}). Thus, we can 
assume that f l=a has affine support — otherwise, we are finished by induction. Similarly, by 
LemmaEl #CSP({/ 4= *}) <t #CSP({/}). Thus we can assume that f l= * has affine support 
— otherwise, we are finished by induction. 

Now, recall that Rf is not affine. Consider any a,b,c £ Rf such that d = a © b © c ^ Rf. 
We have 4 cases. 

Case 1: There are indices 1 < i < j < k such that (aj,6j,Cj) = (aj,bj,Cj): 
Without loss of generality, suppose i = 1 and j = 2. Define the function /' of arity (k — 1) by 
f'(r2, • • • , ffc) = f(r2, ^2, ■ ■ ■ , Tk)- Note that Rf is not affine since the condition o©6ffic ^ Rf is 
inherited by Rf. So, by induction, #CSP({/'}) is #P-hard. Now note that #CSP({/'}) < T 
#CSP({/}). To see this, note that any instance I\ of #CSP({/'}) can be turned into an 
instance I of #CSP({/}) by repeating the first variable in the sequence of variables for each 
constraint. 

Case 2: There is an index 1 < % < k such that a% = bi = q: Since d is not in Rf and 
di = a,i, we find that f l=ai does not have affine support, contrary to earlier assumptions. 

Having finished Cases 1 and 2, we may assume without loss of generality that we are in 
Case 3 or Case 4 below, where {a, (3} £ {0, 1}, a = 1 — a, (3 = 1 — f3 and a', b' , d £ {0, l} fc ~ 2 . 

Case 3: a = (at, (3, a'), b = (a,f3,b ! ), c = (a,(3,c'): Since Rf=* is affine and a, b 
and c are in Rf, we must have either d = (a,(3,d') £ Rf or e = (a,f3,d') £ Rf, where 
d! = a! ®b' © d ' . In the first case, we are done (we have contradicted the assumption that 
d Rf), so assume that e £ Rf but d Rf. Similarly, since Rp=* is affine, we may 
assume that g = (a,f3,d') £ Rf. Since Rp=a is affine and a, b and e are in Rf, we find 
that h = a © b © e = (a,(3,d) £ i?/. Since Rp=$ is affine and a, c and g are in Rf, we 
find that i = (a,(3,b') £ i?/. Also, since Rp = p is affine and a, h and i are in Rf, we find 
that j = (a,/3,d') £ ify. Let f'(r\,r2) = f(ri, r-i, ds, . . . , d^). Since e, g and j are in ify 
but d is not, we have (at, (3), (a, (3), (a, (3) £ Rf, but (a,/5) ^ i?/'. Thus, /' does not have 
affine support and #CSP({/'}) is #P-hard by induction. Also, #CSP({/'}) < T #CSP({/}) 
by Lemma O 

Case 4: a = (a, a, a'), b = (a,a,b r ), c = (a,a,d): Since is affine and a, b and c 

are in Rf but d is not, we have e = (a, a, d') £ Rf. Similarly, since Rf2=* is affine and a, b 
and c are in Rf but d is not, we have g = (a,a,d') £ Rf. Now since Rji=a is affine and a, b 
and e are in ify, we have h = (a, a, d) £ Rf. Also, since Rp= a is affine and a, b and g are in 



10 



Rf, we have i = (a, a, d) G Rf. 

Let f'[T\,T2) = f(ri, r2, C3, . . . , Cfc). If j = (a,a,c') i?/ then /' does not have affine 
support (since c, h and i are in Rf) so we finish by induction as in Case 3. Suppose j G ify. 
Since Rti=& is affine and a, 6 and j are in Rf, we have £ = (a, a, df) G i?/. Let /"(n, ^2) = 
/(ri, r 2 , ^3, . . . , dfc). Then /" does not have affine support (since e, g and £ are in Rf but d 
is not) so we finish by induction as in Case 3. □ 

Lemma [TT] showed that #CSP({/}) is #P-hard when / does not have affine support. The 
following lemma gives another (rather technical, but useful) condition which implies that 
#CSP({/}) is ^tP-hard. We start with some notation. Let / be an arity-fc function. For a 
value b G {0, 1}, an index i G {1, . . . , k}, and a tuple y G {0, l} k ~ l , let y l=b denote the tuple 
x G {0, l} k formed by setting Xi = b and Xj = yj 

We say that index i of / is useful if there is a tuple y such that f(y l ~°) > and f{y l ~ l ) > 0. 
We say that / is product-like if, for every useful index i, there is a rational number Aj such 
that, for all y G {0,1}* _1 , 

f(y l=0 ) = \f{y t=l ). (2) 

If every position i of / is useful then being product-like is the same as being of product type. 
However, being product-like is less demanding because it does not restrict indices that are 
not useful. 

Lemma 12. If f G is not product-like then #CSP({/}) is #P-hard. 

Proof. We'll use Corollary [9] to prove hardness, following an argument from [9]. Choose a 
useful index i so that there is no Aj satisfying ([2]). 

Suppose / has arity k. Let A be the 2 x 2 k ~ 1 matrix such that for b G {0, 1} and 
y G {0, A b)V = f(y i=b ). Let A 1 = AA T . 

First, we show that Eval(^4') is #P-hard. Note that A' is the following symmetric 2x2 
matrix with non- negative rational entries. 

( EyAl y EyA ,yA 1 , y \ = ( Eyfif^f Ey f (V^) f (V^) \ 

U s Vi, 9 EyA 2 hy J \Eyf(y l=0 )f(y l=1 ) E y f(y l=1 ) 2 J 

Since index % is useful, all four entries of A' are positive. To show that 'Eval(A') is #P-hard 
by Corollary [H we just need to show that its determinant is non-zero. By Cauchy-Schwartz, 
the determinant is non-negative, and is zero only if Aj exists, which have assumed not to be 
the case. Thus Eval(^4') is #P-hard by Corollary 

Now we reduce Eval(A') to #CSP({/}). To do this, take an undirected graph G which 
is an instance of Eval(^4'). Construct an instance Y of #CSP({/}). For every vertex v of G 
we introduce a variable x v of Y. Also, for every edge e of G we introduce k — 1 variables 
x e! i, . . . , £ e fc_i of Y. We introduce constraints in Y as follows. For each edge e = (v, v') of G 
we introduce constraints f(x v ,x ei x, . . . ,x ei k-i) and f(x v >,x ei i, . . . ,x e ^_i) into Y, where we 
have assumed, without loss of generality, that the first index is useful. 

It is clear that EvAh(A') is exactly equal to the partition function of the #CSP({/}) 
instance Y. □ 

For w G Q + , let U w denote the unary function mapping to 1 and 1 to w. Note that 
Uq = So, and U± gives the constant (0-ary function) 1, occurrences of which leave the partition 
function unchanged. So, by Lemma El we can discard these constraints since they do not 
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add to the complexity of the problem. Note, by the observation above about proportional 
functions, that the functions U w include all unary functions except for S± and the constant 0. 
We can discard 5± by Lemma [51 and if the constant function is in J-, any instance / where 
it appears as a constraint has Z(I) = 0. So again we can discard these constraints since they 
not add to the complexity of the problem. 

Thus U w will be called nontrivial if w ^ {0, 1}. Let ©^ : {0, l} k — ► {0, 1} be the arity-A 
parity function that is 1 iff its argument has an odd number of Is. Let ->(Bk '■ {0, l} fc — > {0, 1} 
be the function 1 — The following lemma shows that even a simple function like ©3 can 
lead to intractable #CSP instances when it is combined with a nontrivial weight function U\. 

Lemma 13. #CSP(©3, U\, So, 5\) and #CSP(->©3, U\, So, 81) are both #P -hard, for any pos- 
itive A 7^ 1. 

Proof. We give a reduction from computing the Weight Enumerator of a Linear Code, which 
was shown to be #P-hard in Lemma PTOl In what follows, it is sometimes convenient to view 
fflfc, So, etc., as relations as well as functions to {0, 1}. 

We first argue that for any k, the relation ©^ can be simulated by {©3, <5o, <5i}. For 
example, to simulate x\ © • • • ©:% for k > 3, take new variables y, z and w and let m = [A/2] 
and use x\ © • • • © x m © y and x m+ i © • ■ ■ © © z and y © z © w and So(w). 

Since {©3,£o>^i} can be used to simulate any relation ©&, we can use {©3,5o,5i} to 
simulate an arbitrary system of linear equations over GF(2). In particular we can use them 
to simulate the subspace T of code words for a given generating matrix A. 

Finally, we can use U\ to simulate the function which evaluates the weight enumerator on 
T. Then, since A ^ 0, 1, we can apply Lemma[10]to complete the argument. The same proof, 
with minor modifications, applies to -i©3. □ 

Lemma 14. Suppose /E^is not of product type. Then, for any positive A 7^ 1, there exists 
a constant c, depending on f , such that #CSP({f,6o,Si,U\,U c }) is #P-hard. 

Proof. If / does not have affine support, the result follows by Lemma [TT1 So suppose / has 
affine support. Consider the underlying relation Rf, viewed as a table. The rows of the table 
represent the tuples of the relation. Let J be the set of columns on which the relation is not 
constant. That is, if i 6 J then there is a row x with x% = and a row y with yi = 1. Group 
the columns in J into equivalence classes: two columns are equivalent iff they are equal or 
complementary. Let k be the number of equivalence classes. Take one column from each of 
the k equivalence classes as a representative, and focus on the arity-A; relation R induced by 
those columns. 

Case 1: Suppose R is the complete relation of arity k. 

Let /* be the projection of / onto the k columns of R. By Lemma [6l 

#CSP({/*}) < T #CSP({/}) < T #CSP({/,5oA,t/ A ,[/ c })- 

We will argue that #CSP({/*}) is ^P-hard. To see this, note that every column of /* is 
useful. Thus, if /* were product-like, we could conclude that /* was of product type. But 
this would imply that / is of product type, which is not the case by assumption. So /* is not 
product-like and hardness follows from Lemma [T2l 

Case 2: Suppose R is not the complete relation of arity k. 

We had assumed that Rf is affine. This means that given three vectors, x, y and z, in Rf, 
x © y © z is in Rf as well. The arity-A relation R inherits this property, so is also affine. 
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Choose a minimal set of columns of R that do not induce the complete relation. This exists 
by assumption. Suppose there are j columns in this minimal set. Observe that j ^ 1 because 
there are no constant columns in J. Also j ^ 2, since otherwise the two columns would be 
related by equality or disequality, contradicting the preprocessing step. The argument here 
is that on two columns, R cannot have exactly three tuples because it is affine, and having 
tuples x, y and z in would require the fourth tuple x © y © z. But if it has two tuples then, 
because there are no constant columns, the only possibilities are either (0,0) and (1,1), or 
(0, 1) and (1,0). Both contradict the preprocessing step, so j > 3. 

Let R' be the restriction of R to the j columns. Now R' of course has fewer than 2 J rows, 
and at least 2 J_1 by minimality. It is affine, and hence must be ©j or -i©j. To see this, first 
note that the size of R' has to be a power of 2 since R' is the solution to a system of linear 
equations. Hence the size of R' must be 2 J . Then, since there are j variables, there can 
only be one defining equation. And, since every subset of j — 1 variables induces a complete 
relation, this single equation must involve all variables. Therefore, the equation is (Bj or —>(Bj- 

Let /' be the projection of / onto the j columns just identified. Let /" be further obtained 
by pinning all but three of the j variables to 0. Pinning j — 3 variables to leaves a single 
equation involving all three remaining variables. Thus Rf» must be ©3 or -i©3. 

Now define the symmetric function /'" by 

f'"(a, b, c) = f"(a, b, c)f"(a, c, b)f(b, a, c)f"(b, c, a)f"(c, a, b)f'(c, b, a), 

Note that Rf" is ©3 or ->03, since Rf is symmetric and hence Rf" = Rp>. 

To summarise: using / and the constant functions So and Si, we have simulated a function 
/"' such that its underlying relation Ryu is either ©3 or -i©3. Furthermore, if triples x and y 
have the same number of Is then f"'{x) = f"'{y). 

We can now simulate an unweighted version of ©3 or -^©3 using /'" and a unary function 
U c , with c set to a conveniently-chosen value. There are two cases. Suppose first that the 
affine support of /'" is -i©3. Then let wq denote the value of f" when applied to the 3-tuple 
(0,0,0) and let w 2 denote /"'(0, 1, 1) = /'"(1,0,1) = / W (1,1,0). Recall that f"'(x) = 
for any other 3-tuple x. Now let c = {wq/w2) 1 ^ 2 ■ Note from the definition of that 
wo and u>2 are squares of rational numbers, so c is also rational. Define a function g of 
arity 3by 5 (a,/3, 7 ) = C/ c (a)C/ c (/3)[/ c ( 7 )/ w (a, /3, 7). Note that g(0, 0, 0) = w and g(0, 1, 1) = 
g(l, 0, 1) = g(l, 1, 0) = c 2 u)2 = Wq. Thus, g is a pure affine function with affine support -i©3 
and range {0, wo}. The other case, in which the affine support of f" is ©3, is similar. 

We have established a reduction from either #CSP(©3, U\, do, S\) or #CSP(-i©3, U\, 5q, Si), 
which are both ^P-hard by Lemma [T3l □ 

Lemma 15. If f G !F% is not of product type, then #CSP({f,So,Si,U\}) is #P -hard for any 
positive A 7^ 1. 

Proof. Take an instance I of #CSP({f,5o,5i,U\,U c }), from Lemma 1144 with n variables 
x\,X2, ■ ■ ■ ,x n . We want to compute the partition function Z(I) using only instances of 
#CSP({/, <5o, Si, U\}). That is, instances which avoid using constraints U c . For each i, let m; 
denote the number of copies of U c that are applied to Xi, and let m = Yl7=i mi - Then we can 
write the partition function as Z{I) = Z(I; c) where 

Z(I;w)= Yj Z{a) \\ w m * = Z(a)w^ miCTi , 

<re{0,l} n i:<Tj=l ffG{0,l} n 
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where Z(a) denotes the value corresponding to the assignment cr(xj) = o~i, ignoring constraints 
applying U c , and w is a variable. So Z(a) is the weight of a, taken over all constraints other 
than those applying U c . Note also that Z(I;w) is a polynomial of degree m in w. We can 
evaluate Z(I; w) at the point w = A J by replacing each U c constraint with j copies of a U\ 
constraint. This evaluation is an instance of #CSP({/, 5q, Si, U\}). So, using m different 
values of j and interpolating, we learn the coefficients of the polynomial Z(I;w). Then we 
can put w = c to evaluate Z(I). □ 

Lemma 16. Suppose f G T% is not of product type, and g 6 J-% is not pure affine. Then 
#CSP({/,<7,5 ,*i}) is#P-hard. 

Proof. If g does not have affine support we are done by Lemma [Til So suppose that g has 
affine support. Since g is not pure affine, the range of g contains at least two non-zero values. 

The high-level idea will be to use pinning and bisection to extract a non-trivial unary 
weight function U\ from g. Then we can reduce from #CSP({/, So, S±, U\}), which we proved 
^P-hard in Lemma [TBI 

Look at the relation R g , viewed as a table. If every column were constant, then g would 
be pure affine, so this is not the case. Select a non-constant column with index h. If there 
are two non-zero values in the range of g amongst the rows of R g that are in column h 
then we derive a new function g' by pinning column h to 0. The new function g' is not 
pure affine, since the two non-zero values prevent this. So we will show inductively that 
#CSP({/,£/',<5o,<5i}) is #P-hard. This will give the result since #CSP({/, g', S , Si}) trivially 
reduces to #CSP({/, g, So, Si}). 

If we don't finish this way, or symmetrically by pinning column h to 1, then we know 
that there are distinct positive values wo and w\ such that, for every row x of R g with in 
column h, g(x) = u>o and, for every row x of R g with 1 in column h, g(x) = w\. Now note 
that, because the underlying relation R g is affine, it has the same number of 0's in column h 
as l's. This is because R g is the solution of a set of linear equations. Adding the equation 
Xh = or Xh = 1 exactly halves the set of solutions in either case. We now project onto the 
index set {h}. We obtain the unary weight function U\, with A = wi/wq, on using the earlier 
observation about proportional functions. This was our goal, and completes the proof. □ 

Lemma [5] now follows from Lemma [8] and Lemma [TBI completing the proof of Theorem [TJ 
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4 Appendix 



The purpose of this appendix is to prove Lemma [7J for an arbitrary fixed domain [q] . We 
used only the special case q = 2, which we stated and proved as Lemma However, pinning 
appears to be a useful technique for studying the complexity of #CSP, so we give a proof of 
the general Lemma [3 which we believe will be applicable elsewhere. 

Lemma E For every T 'C T q , #CSP(T U |J ceM Q <t #CSP(J r ). 

In order to prove the lemma, we introduce a useful, but less natural, variant of #CSP. 
Suppose J 7 C T q . An instance 7 of #CSP^(J 7 ) consists of a set V of variables and a set C of 
constraints, just like an instance of #CSP(J ?r ). In addition, the instance may contain a single 
extra constraint C applying the arity-g disequality relation \+ with scope (i>c,i> • • • i v C,q)- 

The disequality relation \^ is defined by x^( x ii ■ • • > x q) = 1 if x x, . . . , x q G [q] are pairwise 
distinct. That is, if they are a permutation of the domain [q\. Otherwise, x^( x i> ■ ■ ■ > x g) = 0- 

Lemma [7] follows immediately from Lemma IT71 and [181 below. 

Lemma 17. For every T C JF g , #CSP(^"U lJ ceM 4) <T #CSP^(.T). 

Proof. We follow the proof lines of Lemma EJ but instead of subtracting the contribution 
corresponding to configurations in which some tj's get the same value, we use the disequality 
relation to restrict the partition function to configurations in which they get distinct values. 

Say that J- is symmetric if it is the case that for every arity-fc function / G T and every 
tuple x G [q] k and every permutation it : [q] — > [q], f(x\, . . . , Xk) = f(n(xi), . . . , ir{x k )). 

Let 7 be an instance of #CSP(.F U U c e[g] ^c) with variable set V. Let V c be the set of 
variables u € V to which the constraint 5 c (v) is applied. Assume without loss of generality 
that the sets V c are pairwise disjoint. Let V q = V \ U c e[<j] ^c- We construct an instance I' of 
#CSP^(^ r ). The instance has variables V q U {to, . . . , Every constraint C of 7 involving 

a function f € J- corresponds to a constraint C" of V . Here C" is the same as C except that 
variables in V c are replaced with t c , for each c G [q]. Also, we add a new disequality constraint 
to the new variables to, ... , t q -\. 
Case 1. J- is symmetric: 

By construction, Z(I') = Y^ ,-,Vi-i ' = y °' ' ' ' ' <7 (^- 1 ) = where the sum 

is over all permutations yo, . . . , y q -\ of [5]. By symmetry, the summands are all the same, so 
Z(I') = q\Z{V I cx(io) = 0, . . . , u{t q - x ) = q - 1) = g!Z(i). 
Case 2. is not symmetric: 

Say that two permutations tt\ : [g] — » [q] and 7r2 : [g] — > [g] are equivalent if, for every 
/ G f and every tuple a; G [g] fc , /(-7Ti(xi), . . . , 7Ti(a;&)) = /(^(iCi), . . . , ^(ajfc)). Partition the 
permutations it : [g] — > [g] into equivalence classes. Let /1 be the number of equivalence classes 
and rii be the size of the i'th equivalence class, so n\ + • ■ ■ + = q\^ Let {tt\, . . . , 71^ } be 
a set of representatives of the equivalence classes with tt\ being the identity. We know that 
n\ 7^ q\ since T is not symmetric. 

For a positive integer t we will now build an instance Ii by adding new constraints to I'. 
For each tti other than tt\ we add constraints as follows. Choose a function fi G T and a tuple y 
such that fi(yi, . .. ,y k ) / fi(iTi(yi), . . .,iTi{y k )). If ...,y k )> /^(yi), . . . ,TTi(y k )) then 

define the fc-tuple x l by (a^, . . . , xji.) = {yi, ■ ■ ■ ,y k )- Otherwise, let n be the order of the 

6 In fact, it can be shown that these equivalence classes are cosets of the symmetry group of /, and hence 
are of equal size, though we do not use this fact here. 
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permutation 7Tj and let g r denote /i(7if(yi), ■ ■ ■ ,7r[(yfc)). Since go < g\ and g n = go there 
exists a £ G {1, . . . , n — 1} such that g^ > g^+i- Let (x\, . . . , x^) = (tp (yi), . . . , ^(yk)) so 

• • • ,4) > /i(f*0«a)v • • '^(4)); 
Let u)jj denote fi(iTj(x\), . . . , 7Tj(xJ.)) so, since 7Ti is the identity, we have just ensured that 

w%\ > wu. Let s l = (t i ,...,t .), and let < Zi < h (i = 2, . . . , h) be positive integers, which 

i fc 

we will determine below. Add tzi new constraints to I[ with relation fi and scope s l . Let 
Xi = n^=2 u; 7i- Note that, given er(to) = 7r«(0), . . . ,a(t q -i) = -Ki{q — 1), the contribution to 
Z(J'^) for the new constraints is 

h 

7=2 

So 

h 

Z(l[) = m Z{ I' \ <j(to) = 7Ti(0), . . . , = 7n(q - 1) ) Af . 

i=l 

We have ensured that Ai > 0, since wu > wu > 0, so wu > for all i = 2, . . . , h. We now 
choose the z%s so that \ ^ Ai for all i = 2, . . . , h. If = for any 7 = 2, . . . , h, we have 
Ai = and hence Aj 7^ Aj. Thus we will assume, without loss of generality, that w-yi > for 
all 7 = 2, . . . , h and i = 2, . . . ,h' , where h' < h. Then we have 

e £$=3«7<*T (j = 2,...,/i'), 

where a 7 j = ln(t(; 7 j/w 7 i). Note that an < 0, since -w^ < wu. We need to find an integer 
vector z = (Z2, ■ ■ ■ , z^) so that none of the linear forms Ci(z) = ^7=2 ar yi z i ls zero i f° r 
i = 2,...,h'. We do this using a proof method similar to the Schwartz-Zippel Lemma. 
(See, for example, |21j.) None of the Ci(z) is identically zero, since an 7^ 0. Consider the 
integer vectors z E [/i]' 1 " 1 . At most h h ~ 2 of these can make Ci(z) zero for any i, since the 
equation £i(z) = makes Zi a linear function of z 7 (7 7^ i). Therefore there are at most 
(h! — l)h h ~ 2 < h^ 1 such z which make any Ci(z) zero. Therefore there must be a vector 
z G for which none of the jCi(z) is zero, and this is the vector we require. 

Now, by combining terms with equal Aj and ignoring terms with Aj = 0, we can view 
Z(I' e ) as a sum Z(I' e ) = ^ CjAf where the Aj's are positive and pairwise distinct and 

ci = mZ(l' I er(io) = 0, . . . , cr(tg_i) = g - 1). 
Thus, by Lemma 3.2 of [8j we can interpolate to recover c\. Dividing by ni, we get 

Z(I' I ff (t ) = 0, . . . , = 9 " 1) = □ 

Lemma 18. For every ^ C T q , #CSP^(.F) < T #CSP(J r ). 

Proof. We use Mobius inversion for posets, following the lines of the proof of [2J Theorem 8] 
Consider the set of partitions of [q]. Let denote the partition with q singleton classes. 
Consider the partial order in which 77 < iff every class of r/ is a subset of some class 

7 Lovasz |15| had previously used Mobius inversion in a similar context. 
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of 6. Define /j(0) = 1 and for any 6^0 define fx{9) = — J2 v <e n±e M 7 ?)- Consider the sum 
J2 V <8 / i ( r /)- Clearly, this sum is 1 if 6 = 0. From the definition of /j, it is also easy to see that 
the sum is otherwise, since 

J>fa)=M(0)+ E Mfa)=0. 

r)<8 V<S,r]^e 

Now let / be an instance of #CSP^(J r ) with a disequality constraint applied to variables 
to, . . . ,t q —i. Let V be the set of variables of /. Given a configuration a : V — > [g], let i?(<r) 
be the partition of [g] induced by of (a (to), . . . , a(t q -\)). Thus i and j in [q] are in the same 
class of i?(<j) iff <r(ij) = cr(tj). We say that a partition 77 is consistent with a (written 77 =<! a) 
if 77 < ??(cr). Note that 77 ^ a means that for any i and j in the same class of 77, cr(tj) = c(ij). 

Let be the set of configurations a that satisfy all constraints in / except possibly the 
disequality constraint. Then Z(I) = J2aen w ( a )^-<T-> where l a = 1 if a respects the disequality 
constraint, meaning that #(cr) = 0, and 1 CT = otherwise. By the Mobius inversion formula 
derived above, 

z(i) = e^) E tii)- 

aen r;<i?(cr) 

Changing the order of summation, we get 

z(i) = em^E E = E^) E 

Now note that Ylcr-rj^cr w ( a ) is the partition function Z(I V ) of an instance I v of #CSP(J ?r ). 
The instance 1^ is formed from I by ignoring the disequality constraint, and identifying 
variables in to,. . . ,t q -i whose indices are in the same class of 77. Thus we can compute all 
the Z{I n ) in #CSP(.F). Finally, Z{T) = fi(rj)Z(I v ), completing the reduction. □ 
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